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Abstract 



The spectrum of anomalous dimensions of gauge-invariant operators in maximally super- 
symmetric Yang-Mills theory is believed to be described by a long-range integrable spin 
chain model. We focus in this study on its si (2) subsector spanned by the twist-two single- 
trace Wilson operators, which are shared by all gauge theories, supersymmetric or not. 
We develop a formalism for the solution of the perturbative multiloop Baxter equation 
encoding their anomalous dimensions, using Wilson polynomials as basis functions and 
Mellin transform technique. These considerations yield compact results which allow ana- 
lytical calculations of multiloop anomalous dimensions bypassing the use of the principle 
of maximal transcendentality. As an application of our method we analytically confirm the 
known four- loop result. We also determine the dressing part of the five-loop anomalous 
dimensions. 



1 Introduction 



The success of gauge theories in accurately describing the laws of nature is based on the 
availability of computational techniques, see e.g., Ref. pQ, which allow for a systematic 
improvement of approximations involved. Perturbative expansions in the gauge coupling 
constant g YM are conventionally deduced from Feynman diagrams. However, due to uncon- 
trollable proliferation of the latter at higher orders in g YM , the rules quickly become un- 
manageable, making direct computations already at four-loop order highly nontrivial and 
require massive computer manipulations. On top of this, individual Feynman diagrams 
obscure underlying properties of the theory and reveal simple results enjoying sometimes 
enhanced symmetries only in their sum. One was therefore compelled to search for an 
alternative approach which presented itself recently. 

On the one hand, some time ago it was established that at weak coupling one-loop 
spectra of anomalous dimensions of maximal-helicity gauge-invariant operators in QCD 
coincide with energy spectra of a one-dimensional non-compact Heisenberg magnet 0, [3] . 
The latter can be diagonalized by means of the traditional Bethe ansatz formalism of 
integrable systems and yields anomalous dimensions of the corresponding four-dimensional 
gauge theory. These simplifications are echoed by higher loop contributions, especially in 
super symmetric gauge theories. It was found in Refs. [U [5] that all single-trace operators 
in planar, maximally supersymmetric gauge theory 

O = tr (X(D 2 + X)YZXXXF +± (D + Y)X . . . ) (1.1) 

can be described by a long-range integrable spin-chain model with elementary excitations 
identified with the particle fields Y, Z, A etc. of the gauge theory and/or covariant deriva- 
tives acting on them propagating on the vacuum state |0) = tr (X L y Less supersym- 
metric Yang-Mills theories entertain integrability only in certain closed subsectors under 
renormalization group evolution [BJ. 

On the other hand, the AdS/CFT correspondence [7] conjectures that the strongly 
coupled M = 4 Yang-Mills theory is dual to a free type IIB super-string theory on an 
AdSsxS 5 background. The latter was found to be classically integrable as well jH]. Using 
this conjecture as a virtue led to a suggestion of an integrable structure which interpolates 
between weak and strong coupling regimes. Though the underlying spin chain model is not 
known, a set of Bethe ansatz equations is nevertheless available [5J E], which has passed 
a number of non-trivial tests at weak coupling, see e.g. j9] and [TO] , as well as at strong 
coupling by positive comparison with perturbative string theory, see e.g. [IT] and [12] . 

These findings suggest to use the putative integrable structure as an alternative to 
the conventional Feynman diagrams technique for multiloop calculations of anomalous 
dimensions. In this paper, we develop a practitioner's formalism building up on earlier 
considerations based on the all-order Baxter equation [T3l [T4"] for finding the spectrum of 
twist-two Wilson operators 

O = tr (XD%X) . (1.2) 

These arise in all gauge theories albeit with a different field content, the scalar X being spe- 
cific to supersymmetric cousins of QCD. Their anomalous dimensions have been obtained 
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diagrammatically to a considerably high-order [T5], HH [T7J HH [19] . 

The Baxter equation is advantageous over the Bethe ansatz formalism if one is interested 
in a systematic analytical framework. However they both enter on equal footings for 
numerical studies, and Bethe equations were used in the past together with the principle 
of maximal transcendentality [20] to perform phenomenal computations [211 [22] . 

Our following consideration is a generalization of the study in Ref. [23] which was based 
on a deformation of the solution to the one-loop Baxter equation. What will differ in the 
current work is that we will introduce a new basis of functions used in the construction of 
next-to-leading order solutions, the so-called Wilson polynomials. For comparison we also 
present the basis of continuous Hahn polynomials used in [23]. Furthermore we obtain a 
new form for non-polynomial contributions which is free from multiple sums involving Stir- 
ling numbers. The latter property is essential for obtaining analytical results for anomalous 
dimensions in terms of nested harmonic sums. Our subsequent presentation is organized 
as follows. In the next section, we briefly review the formalism of the Baxter equation in 
maximally supersymmetric gauge theory and then present a novel form of the solution in 
two- to four-loop order in the gauge coupling. The non-polynomial parts of the Baxter 
equation are analyzed in Section [5] using the Mellin transform technique. We present the 
analytic form of anomalous dimensions and then in Section [71 we discuss the reciprocity 
properties including the dressing part of the five-loop anomalous dimensions. Finally, we 
conclude. Several appendices summarize basic definitions required in the main body of the 
paper and details of calculations which are two lengthy to be presented in the main text. 



2 Baxter equation 

The spin-chain description allows one to calculate anomalous dimensions of Wilson oper- 
ators as a function of the 't Hooft coupling constant]]] g 2 = g 2 M N c / (An 2 ). However, the 
formalism based on the Bethe ansatz equations has the drawback that its predictions fail 
when the order of the perturbative expansion in g 2 exceeds the length L of the operator 
under study [2JJ[25]. This implies that for twist-two operators (II. 2p the onset of wrapping 
effects occurs starting from four-loops already and the complete anomalous dimension is a 
sum of two terms 

7(^)=7 (asy) (^) + 7 (wrap) (^)- (2-1) 

The first contribution on the right-hand side is determined by the solution to the 

asymptotic Bethe ansatz equations and can be written in terms of the Baxter function 
Q(u) as [13] 

('"IB)' • 

1 Please note that we use a different convention for the coupling constant then in [23j . 
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The latter is a degree-M polynomial in the spectral parameter u with zeros determined by 
the Bethe roots 



M 



Q(u) = H(u-u k (g)). (2.3) 



k=l 



It obeys an equation known as the asymptotic Baxter equation 

(x + ) L e a+{u+) - e{u+) Q(u + + (x~) V- (u-) - e(u-) Q(« - i) = t{u)Q{u) . (2.4) 

Compared to the Baxter equation for the familiar non-compact nearest-neighbor XXX 
Heisenberg spin chain 

<B[Q ] = (u + ) L Q (u + + (u~) L Q (u - i) - t {u)Q {u) = , (2.5) 

with factors (u^) 1 = (u ± |) L accompanying the corresponding Baxter polynomials, Eq. 
(12. 4p possesses highly non-trivial dressings factors reflecting coupling-constant dependent 
dynamics of the four-dimensional Yang- Mills theory. First, the spectral parameter gets 
renormalized [26] and reads x = x[u] = |(u + a/m 2 — g 2 ), with the assumed conventional 
notation x ± = and, second, the exponents a and provide the interpolation between 

weak and strong-coupling expansions [9] and read 



r=2 s=r+l 71 V 




s ' 2 u s ^(t) 



x r ~ l 




(2.7) 



with the expansion coefficients given by 



0,fa)=tfn(i(.-r))/ ■ (2.8) 

Since the Bethe roots acquire dependence on the 't Hooft coupling, the Baxter function 
can be expanded in a perturbative series Q(u) = Qq(u) + g 2 Qi{u) + . . . and each term 
found explicitly as a solution to Eq. (12.41) as we demonstrate next. Notice that all sub lead- 
ing Baxter functions Qe>o{u) are polynomials in the spectral parameter of a degree two 
units lower than the leading Qq{u). To four- loop order, the dressing functions admit the 
expansion 



° ± (u) = g 2t -it + 9 4 



(2.9) 



2 This nomenclature should not be confused with the dressing factor related to the phase Q(u), which 
we also refer to in later sections. 
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jit - ± (Hit + HP - - i(7 ± ) 2 « ± + x*) 

Ub ill 

+ . 



e(w) = (/ 6 C3ate 



+ ..., (2.10) 



t{u) 



where the expansion coefficients are introduced explicitly in Appendix lAl 

Yet another unknown in Eq. (12.41) is the transfer matrix, which takes the form 

Ke(x+) L (2 + J2&k(g)Mx + r k ) -^%( 9 )3m(x + )^. (2.11) 

V fc>l / k>l 

Here the upper limits in the sums can exceed the length of the operator in question and the 
emerging charges Q.k>L along with D^fc serve to compensate non-polynomial terms arising 
in the left-hand side of the finite difference equation (12. 4ft stemming from the expansion of 
the renormalized rapidity parameter and dressing factors in Taylor series in the 't Hooft 
coupling. The charges admit perturbative expansions 

n fc =4 01 + /4 1] + s 4 4 2] + ..., (2.i2) 
^=4 01 + ^i 1] + ^i 21 + ---- 

And the only non-trivial contributions for L = 2 operators up to four-loop order read 

= -M(M + 1) , (2.13) 
n W = _ (2 M+l)3?e[ 7o + ] , 
Ojjfl = _(2M + l)He [ 7l + ] - |»e [3( 7o +) 2 + «+] , 

£jM = _(2Af + l)«e [ 7 +] + |«e [a + (( 7o + ) 2 + « + ) - 7o + (67i + + 8/3+ + C 3 ) - 4 X + ] , 

J3] 



^ J = -^e[ 7o + (( 7o + ) 2 + « + )] 

[3] 
2 

[3] 



9t£ J = 



= £8te [To + ] • 

Here 0^ * s the eigenvalue of the quadratic Casimir operator of the collinear conformal 
subgroup in the basis of conformal Wilson operators such that the leading order transfer 
matrix admits the conventional form for the two site non-compact Heisenberg spin chain, 

t ( u ) = ( u + f + (u-f + 2 01 . (2.14) 

Finally, the solution to the Baxter equation has to be supplemented with the condition of 
the vanishing quasi-momentum 

w= ir * h «±h*> o, (2 . 15 ) 



7T 



iv^ Q(-l-gt) 



in order to pick out only cyclic, physical states. 
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3 Wilson vs. Hahn 



It is known for quite some time that the leading order solution Qq for the non-compact 
two-site Heisenberg magnet is given by the continuous Hahn polynomials^ [28] . 



Qo{u) 



3-^2 



-M, M + 1, ± + iu 
1, 1 



(3.1) 



Proceeding to higher loops, it was demonstrated in Ref. [23] that subleading contributions 
Qe>o{ u ) to the Baxter function can be obtained by a deformation of the leading order 
result (13.11) . The equations which these corrections obey remain of second order in finite 
differences, but acquire inhomogeneous terms depending on lower-order functions. This 
implies that the structure of all polynomial higher-loop contributions can be immediately 
understood once the building blocks for the two-loop Baxter function are known. To find 
the latter it suffices to expand both side of the Baxter equation (12.41) to 0(g 2 ) and find 



= [(« 



1„2 
2 



g 2 (l + il0 u')]Q(u-t) (3.2) 
-<? 2 + jQ 2 ] Q{u), 



h 2 lo)(M 



ig 2 7o) renormalized 



with the quadratic conformal Casimir H 2 ~ — (M + 1 
by the one-loop anomalous dimension 7 = 2S\(M) to this order of perturbation theory. 
Matching this to the equation obeyed by the continuous Hahn polynomials, we obtain the 
result for the Baxter function with incorporated two-loop corrections 



Q{u) = N(g) 3 F 2 



-M, M + 1 + 3 2 7o , \ + iu + j-g+ ^ 2 7o 
1 + iy/2g + i(7 2 7o , 1 + ^ 2 7 



1 



(3.3) 



Expanding this in a Taylor series with respect to the 't Hooft coupling, we find for the 
two-loop correction itself 



Qi(u) = hQ (u) + \ (2 l0 d Sl - d\ - dl) 3 F 2 



-M, M + l + 2<Ji, ~+iu + 8 2 
l + *i + *2 + 53, l + <5 2 -<5 3 



<5,=0 

(3.4) 

where we chose a spin-dependent form of the perturbative expansion of the normalization 
constant N(g) ~ 1 + g 2 bi(M) with 



6i(M) ^AS 2 + S 2 -2S 1 S 1 



(3.5) 



in order to reduce the degree of Q\ in accordance with the definition (12. 3p such that 
deg Qi = (deg Qq — 2). Here and below the nested harmonic sums (1A.4I) appear as functions 
of two arguments 

Sai,a,2,... — *S'ai,a2,...( J ^) ) *S'ai,a2,... = | S'oi,02,... (2Af) . 



3 We summarize their basic properties in Appendix [Bj 
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The Baxter function for twist-two operators in the basis of continuous Hahn polynomials 
has been obtained analytically to three- loop order in [23]. However, as can be seen from Eq. 
( 13 .41) the number of deformation terms arising is quite substantial. This calls for a quest to 
find a more concise representation. To this end one notices that the Baxter function for the 
ground state is symmetric under the reflection u — > —u, which is however not transparent 
in the representation in terms of Hahn polynomials (13. ip but becomes explicit in the basis 
of the Wilson polynomials^ 



Qo(u) 



4-^3 



M M+l 
2 ' 2 



, 2 + ill, 2 

1 I 

" ' ' 2 



IU 



(3.6) 



Analogously to the previous consideration, matching the two-loop Baxter equation to the 
equation for Wilson polynomials (IB. 7|) . we find the two-loop solution 



Q(u) = N(g) 4 F 3 



M 
' 2 ' 



^ + !<7 2 7o, § + ««, \ 



IU 



(3.7) 



up to an overall coupling-dependent normalization constant N(g). Expanding this result 
to order g 2 yields the two-loop Baxter polynomial Q\{u) in Wilson basis 



Qi(u) = a 1 {M)Q {u) + \{2 l0 d Sl - <9|) 4 F 3 



M M+l 
2 ' 2 



+ 5 1 



' 2 



1 + iu, h — iu 



l + di + 52, 1 - 52, 

where the normalization constant is N(g) ~ 1 + g 2 a\{M) with 

ai (M) = 35? + S 2 + 5_ 2 - 2S 1 Si . 



2 
1 
2 



$1,52=0 

(3.5 



(3.9) 



A few comments are in order concerning the relation of this representation to the one in 
terms of continuous Hahn polynomials. First of all, there is one deformation less. As 
a consequence, the number of polynomial contributions in higher loops will drastically 
decrease. Second, the deformed parts in the two representations are not identical and 
hence their degree reduction coefficients multiplying the leading order solution differ as 
well. However, it should be noted that these degree reduction coefficients can, in both 
cases, not contribute to the anomalous dimension by symmetry arguments. Although less 
obvious, the same is true for higher-loop contributions. 

Let us now turn to higher loops contributions and use the finding of this section to 
devise an efficient formalism to determine perturbative solutions. We will present the 
results both in the basis of Wilson and continuous Hahn polynomials. 



4 Polynomial contributions 

To start with, as we observed in the previous section in order to tackle higher order 
corrections to the Baxter function it suffices to introduce a doubly-deformed function and 

4 We briefly review them in Appendix [Bj 
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its derivatives with respect to the deformation parameters 



M M+l , s 1 _i_ An. 1 —im 

5l 52 3 1 l + 5x + 5 2 , 1-S 2 , I 



(4.1) 

<5 1= <5 2 =0 



Then a straightforward scheme presents itself for the construction of the £— th order func- 
tion Qi(u). That is, Q^(it) is a linear superposition of the structures T^ ^) with £x + £ 2 < ^ 
accompanied by degree-<i transcendental numbers of d = 2£—2£ 2 —£i and a QcrP ro P or ti ona l 
term, such that the degree of Qi is reduced to degQ^ = (degQ — 2). 

Analogously, in the basis of continuous Hahn polynomials we have a triple- deformed 
function 



-M, M + 1 + 5 U \ + iu + 5 2 
1 + S x + 5 2 + 5 3 , l + 5 2 -5 3 



■ (4-2) 

<5i=i52=53=0 



In the basis of Hahn polynomials there are more possible deformations than in the 
Wilson basis. As a consequence the number of terms at each order of the perturbative series 
is also considerably increased. As was shown in Ref. [23] at three loops, a simple counting^) 
gives a total number of eleven terms. Changing the representation from 3 F 2 —>■ 4F3 reduces 
the number of contributions to the Baxter function by almost half. At four-loop order 
this effect will decrease the number of contributions from 36 to 15 terms, as demonstrated 
below. 

In the following we will divide the contributions to the Baxter functions into polynomial 
and non-polynomial contributions 

Q e (u) = Q ( f\u) + Q i ; p \u) } (4.3) 

respectively. The terminology used here needs clarification. Of course, at any given order of 
perturbation theory, the Baxter equation is polynomial. However, it consists of two types 
of terms, the first one explicitly polynomial in the spectral parameter u and the other 
containing inverse powers of the spectral parameter m ± accompanying Baxter polynomials 
and thus appearing superficially non-polynomial. Indeed the inverse powers of u conspire 
to cancel in the sum of the latter such that the net result is polynomial as it should. 
However we choose to split the Baxter function according to this nomenclature inherited 
from their source in the equation. 

4.1 Wilson basis 

Following the strategy outlined above we find first the polynomial part of the perturbative 
Baxter function in the Wilson basis. First, the three-loop Baxter function reads 

Q 2 \u) = a 2 Q {u) + ! (71 + CI170) T ( i j0 ) -\(K 2 + 2a x ) T (0 ,i) 



5 Note, that for this representation, there is a term X^o^), which has a third- order derivative w.r.t. 
the deformation parameter. 
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+ §7o T (2.o) - §7o7(i,i) + ^T {0 ,2) , (4.4) 

where the transcendental coefficient K 2 of degree two is a linear superposition of anomalous 
dimensions and inhomogeneities introduced in Appendix |A] 

^2 = « + |7 2 = ho-^-2, (4-5) 

and the normalization function a 2 = a 2 (M) reducing the degree of the higher-loop poly- 
nomial depends on the non-polynomial contribution computed later in Section [5j 

At four-loop order a further transcendental function arises from the Baxter equation, 
see Appendix |AJ from the expansion coefficients of the dressing factors (12. 9p . They appear 
in a certain combination with a degree of transcendentality four, 

K A = -6a 2 + 24 X - e lo + H + | To 4 + 10 7 o7i • (4.6) 
So that finally the polynomial part of the four-loop Baxter function is 

Q ( t ( u ) = a 3 Qo(u) + |(72 + ai7i + ^27o) 7(i,o) - 5 (7Q + 3C37o + 6a x K 2 + 12a 2 ) T (0 ,i) 
+ I (27i7o + ai7o) T (2,o) - f (71 + a i7o + |^27o) ^(1,1) + ik (To + 2K 2 + 2a i) T (o,2) 

+ 4^^0 T (3,0) - 3l7o T (2,l) + 1^270^(1,2) ~ 5^0^(0,3) , (4.7) 

with the degree-reducing coefficient a 3 . In the last equation, the term proportional to £3 
stems from the dressing factor. 

4.2 Hahn basis 

Let us also include for completeness the three- loop Baxter function obtained in [23j, con- 
verted to the notation of Eq. (14. 2p and rescaled coupling, 

Q { 2 p \u) = b 2 (M)Q (u) + |(ti + &i7o) T (1 , ,o) - § (K 2 + 260 (^0,2,0) + T (0 ,o, 2) ) 

+ |7o T (2,0,0) - 87o(T(l,2,0) + 7(1,0,2) + 1^(0,0,3)) 

+ ^(6T( ,2,2) + T(o,4,0) + ^(0,0,4)) • (4.8) 

The appearing functions are the two-loop normalization constant for the Hahn basis &i 
given in (13.51) and the transcendental function K 2 in (14.51) . We omit the precise structure 
of the normalization constant b 2 (M), as it will be of no use for us. 

In the basis of continuous Hahn polynomials the four-loop result is quite lengthy and 
reads in the conventions of Eq. (14.21) . 

= hQo{u) + \{l2 + 6171 + 6 2 7o) 7(i,o,o) 
-i (K A + 3C 3 7o + ^ X K 2 + 1260 P(o A o) + ^0,0,2) ) 

+ 1§2^ 2 (7(0,0,4) + 7(0,4,0) + 6T( 0j2i 2)) 

-^70 [K 2 - 271 + 260 (7 1 (i,2,o) + §7(0,3,0) + 7(i,o,2)) 
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+ 3^0 ( T (l,2,2) + | T (1,0,4) + §2(0,3,2) + | T (1,4,0) + TE T (0,5,0)) 

_ ^7o (2(2,0,2) + 2(2,2,o) + 12(0,4,0) - |2(o,o,4) + §2(1,3,0)) + 3g7o2(3,o,o) • (4.9) 



Again, we omit the definition of the normalization 6 3 . is determined by (14. 6p . These 
expressions for the three- and four-loop Baxter functions ( 14.8[) and ( 14.91) in the Hahn basis 
should be compared to their Wilson basis counterparts in Eqs. (14.41) and (14.71) . respectively. 
It is apparent that the Wilson basis enormously simplifies higher- loop computations and 
should be the method of choice for subsequent loop-orders. 

5 Non-polynomial contributions 

In order to complete the solution to the Baxter equation we have to address the non- 
polynomial contributions. According to the nomenclature of Section HI we collectively 
label all inhomogeneities of the form (u ± )~ k Qg(u) with k > 1 as non-polynomial. Recall 
however, that the transfer matrix ( 12. lip is chosen in such a way, that it compensates these 
non-polynomialities and in the sum of all these contributions at a given loop-order the 
polynomiality is restored. 

Non-polynomial inhomogeneities in the Baxter equation of twist-two operators appear 
for the first time at three-loop order, while wrapping effects set in at four loops. According 
to our choice of splitting in ( 14. 3 p we will complete the three-loop Baxter function by 
obtaining the term Q^S^ * n an nove l form. The non-polynomial contributions at three- 
loop have already been found in [23] • However, the representation used there is given in 
terms of Stirling numbers, which complicates the computation of anomalous dimension. 
Therefore we will present a novel representation of these terms here, which is based solely 
on Mellin transform techniques of the (u ± )~ k with k > 1. The procedure is different to the 
one presented in [23J , as there the effective polynomial of all non-polynomial terms has been 



is that we give the Mellin transform of all single non-polynomialities. This results in a 
toolbox, which allows to simply construct the solution from a general set of expressions. 
The advantage is, that the final representation does not depend on complicated coefficient 
functions involving Stirling numbers, rendering the computation of anomalous dimension 
more feasible. 

Before we turn to the explanation of the method used in our calculation, let us sum- 
marize the results of our analysis in this section. The non-polynomial contributions to the 
three- and four-loop Baxter functions are given by the following expressions 




M 




(5.1) 



k=0 



6 This procedure is given as Lemma 1 in [23] 
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M 



Qf p \u) = Y, 2Ste[ft(«)] R k (M) (r^\M)+r^\M)) , (5.2) 



k=0 

with expansion coefficients determined in Eqs. (15.131) . and (15.391) . (I5.50p . (15.511) respectively. 
These expressions complete the three- and four-loop Baxter function, which is given by 
the sum of (14.41) and ( 15. 1H . and (14.71) and (15. 2p . respectively. The functions used in the 
representation of the non-polynomial parts are given by 

Rk{M) ~ (k\y r(M + i-k)' Pk{u) - k\r$ + iu) ■ (5 - 3) 

Details of the computation are presented in the following subsections, while the complete 
dissection of the Mellin techniques is deferred to Appendix [Dl 

5.1 Three loops 

At three loops order, the non-polynomial contribution 1/2(11) to the Baxter equation 

&[Q 2 ]='-- + U 2 [Qo], (5.4) 

reads 

U 2 [Qo) = [Qo(u + i) - Qo{u)) + ^:Qo(u + i) + cc. . (5.5) 

The procedure for finding a closed solution to this equation via the Mellin transform 
technique together with a complete set of building blocks required for generic higher order 
analysis of non-polynomial contributions is presented in Appendix [Dl We refer to it for a 
complete list of definitions of objects arising in this calculation. Below we merely assemble 
specific terms at three-loop order. 

Extracting the first, constant term from Qq accompanying 70 in Eq. (15.51) and introduc- 
ing a new function Q Q via Eq. (ID. 731) . we can decompose U 2 into a sum of two terms with 
each of them being separately polynomial. Then, making use of Eqs. (ID.90[) . (ID. 931) and 
(ID.95P we can obtain solutions to Eq. (15.41) stemming separately from both combinations 
of non-polynomial contributions. 

The first inhomogeneity in U 2 is of the form of (ID. 901) with L = 0, i.e., 

^ §o(« +<)-—&(«-<)=<£ 2 ^ e[jPp ' l(M)] R P (M) . (5.6) 

p=i 



Then, matching the right-hand side of this equation to Eq. (ID. 181) provides via Eq. ( ID. 251) 



the contribution of this inhomogeneous term to the Baxter function Q 2 (u), which reads 



M 



J2^MPk(u)]R k (M)S 3 (k), (5.7) 



fc=0 
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up to an overall factor |7 = ±17^ accompanying (15.61) in Eq. (15.51) . 

Now turning to the second inhomogeneity in U2, one notices that it is given as a linear 
combination of Eqs. fHX90l) . flLTMl) and fITJMl) with L = 1 and Qo(±f) = 

2&e[P p _ 1 (u)] 



1 ~ 1 00 

__Q ( M + ,) + __Q ( M _ z ) = -iQ' (|) ^ 



p=M+l 



^ 2 8te[P^!( M )] 
p=l r fc=l 

n , n , 1 n < \ n>( M 23?e[P p _ 1 ( M )] 
<5o(«) + 7— TaQo(«)= -«<3o(-|) 2^ 



\ / v I p=M+X c 

-f ^ tW4M), (5.9) 

and (ITX95|) with L = 



p=i ^ fc=i 



1 -J. = - i f; 2 ^(»)l, ,5.10) 

p=i 

respectively. Then, the sought-after combination of terms is 

(Qo{u + Q (w)) + T^p - *) - Qo(«)) + 2*7o ^ - 7^ J 

= f me[P ^ U)] ±R k (M)Z lip ,k-l) 

p=l P k=l 

-«(««) -%H)) E ^^ + 27o fi!M« 

p=M+l ^ p=l ^ 

And since iQ' (±|) = ±70, the contribution of the infinite series ~ Y^=m+i cance ls between 
the last two terms and the result takes a polynomial form 



^ (Qo(u + 1)- Qo(w)) + (Q {u - i) - Q (u)) + 2«7o ^ - -Lj (5.11) 

p=l P ^ k=l > 



Combining this result with the one found earlier in Eq. (15. 7p multiplied by its proper relative 
coefficient in Eq. (15. 5p and identifying the expansion coefficients in the summand with Bk 
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of Eq. (1D.24|) . we can immediately write down the contribution of both inhomogeneous 
terms to Q 2 as 



M 



16 J2 2 ^ e l p k( u )} MM)rffi(M) , (5.12) 

fc=0 

with , introduced in Eq. (ID. 251) . taking the following explicit form 

r%£\M) = £ (9a,i(k) + V 3)1 (k) + 2 lo (W 3 (k) - S 3 (k)f 

= ± (v 3j i(k) + 2 l0 (W 3 (k) - S 3 (k))) , (5.13) 

determined in terms of the following sums 

(5.14) 

R m (M)Z b (r,m-l), (5.15) 

Rm{M) 

Z 6 _i(r,m-1), (5.16) 

m 

R m {M) Z b {r,m) , (5.17) 

with Z b introduced in Eq. (1D.87|) . Restoring the overall coefficient ^ of (15. lip in Eq. (15.51) . 
we recognize (15.121) as Eq. (15.11) . 



k 

E 

i=i 


1 




l a Ri(M) ' 




k 

E 

r=l 


1 


r 

E 

m=l 


r a Rr(M) 


k 

E 

r=l 


1 


r 

E 

m=l 


r a R r (M) 


k 

E 

r=l 


1 


r 

E 

m=l 


r a Rr(M) 



5.2 Four loops 

At four loops the number of non-polynomial contributions is greater, however, due to the 
perturbative iteration, the contribution of inhomogeneities due to the two-loop Baxter 
function Qi has the same form as in Eq. (15. 5p . i.e., U 2 [Qi]- The complete set of non- 
polynomial inhomogeneities then consists of three terms, 

Q3[Q 3 ] = . . . + U 2 [Qi] + U 3 [Q ] + t/i[Q$ np) ] , (5.18) 

where Ui[Q\ ] arises from non-polynomial contributions to the three- loop Baxter function 
Q 2 computed above and has the form 

Ui[Q ( 2 p) ] = {\- Q { 2 np \u + i) + (| - ryo «-) Q { 2 np \u - i) 

-(l + (2M + l) 7o + )Q^ ) ( M ). (5.19) 
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While U3 is another novel non-polynomial function of the leading order Qq, 

UsiQo] = I ^-^i [Qo(u + Qo(u)) + (Qo(u + 1) + Q (u)) (5.20) 

+ + l) ~ 4^ 7o+ ^ + + (7o+)2 ) } + C ' C - ' 

where 

Bf = 10( 7o ± ) 3 + 47? + 67^ - + e ± . (5.21) 
Let us address all three contributions in turn, starting with the latter. 



5.2.1 Inhomogeneity Us 

Following the methodology developed at three-loop order, we split the inhomogeneity U 3 
into a sum of terms whose series representations can be matched into the generic types 
analyzed in Appendix [D] According to results given there, the inhomogeneity in Eq. (I5.20p 
is a linear combination of Eqs. (ID. 901) . (ID.93P with L = 0, 2, 3 and Qo(if) = 0, namely, 

(«+)< (* ( " + !) " * ( "') + (Py * (ll) ) (5 - 22) 

- f 28ie|Pp - l( " )l £ i«M) * - 1) + }z,(p, *>; 

p=l P k=l 

+ £ 23?e| t l(M>1 {HQo(-i)-«(i))^(rt 

p=Af+l ^ 

+1 (Qo(-I) - OS(i)) £i(p) + I (Qo'(-I) " Qo'(|)) } 



multiplied by the factor of ^. The term 



( ^ H 5 , v Q (m + z) + Qo(w)) - (<5o(w - i) + Qo(«)) (5-23) 



p=l ^ k=l 



- E 2 ^ e[p ;~ lN] K^oH) +o / o(i))^(p) - f (qo(-d +gg(i))} 



p=M+l 

accompanied by ^7o~. And the terms 



-\Qo(u + i) - -Q (u -i) = -V 2 ^ e[jPp " l(M)] i? P (M) , (5.24) 

p=l 
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-Qo(u) - ^Q (u) = £ 2 ^ e[i ^ l( " )] £ R k (M) , (5.25) 

with the factors of j^Bf and — t7q" [(7o~) 2 + « + ] > respectively. Further, for the completion 
of the last three equations we need as well Eq. (ID. 951) for L = 2 times §7(j~, and for L = 
entering with ^5+ - | 7 + [( 7 +) 2 + a+] . 

Combining all these expressions and considering contributions proportional to Z 2 (p), 
we find that the infinite series present in separate terms cancel between the two equations 
(I5.22p and Eq. (ID. 951) for L = 2 and their net result is equal to 

.^±^m Z2(p) . ( , 26) 

P =i p 

At the same time, the non-polynomial coefficients accompanying Z\(jp) vanish since the 
leading order Baxter polynomial is an even function of the spectral parameter. Turning to 
the remaining two infinite-series contributions we deduce that they can be resummed into 
a concise expression such that the total inhomogeneity U3 admits the following form with 
clearly separated polynomial terms 

At 



UslQo] = - & E me[Pp - l{u)] { [ IB+ - 8 7o + (( 7o + ) 2 + <*+)] R P {M) (5.27) 
p=i 

P r 

+ ±Z 3 (p, k - 1) - 2 7o + (^i(p, k - 1) - \Z 2 (p, kyj 

+ 8( 7o + a + -8/3+) +2 1 +Z 2 (p) 



k=l 



+8T„ + «T„ + ) 2 + Q + ) 



Here we used the identity 

1 Z L (p, k - 1) = 1 Z L (p, k) + ^ ^Ezfci) (5.28) 

to simplify intermediate results. Even though there are remaining non-polynomial contri- 
butions, (the last line in Eq. ( 15. 27ft ) that do not cancel on their own, they will after we 
add terms stemming from 6^2 [Qi] as will be demonstrated in the next section following the 
same lines of reasoning as in Section I5TT1 upon the replacement Qq(u) — > Qi(u). 

5.2.2 Inhomogeneity U 2 

To start with, we write the two-loop Baxter polynomial as 

M 

Q x {u) = Y,MP P {y)\Rp{M)r x , p {M) , (5.29) 

p=0 
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with (see Ref. [23] ) 

r 1>p (M) = b(M) + 7o (M) (^(p + M) — S 1 (M) - \ S^p)) - S 2 (p) , (5.30) 

and 

b(M) = foi(M) - S?(M) , (5.31) 

where b\ was introduced in Eq. (I3.5p . such that fi i0 (M) = Then we split Qi into a 

constant piece and the rest Q\{u) starting from P p>0 (u) 

Q 1 (u) = n fi (M)+Q 1 (u). (5.32) 

Then as in Section |5TT| the inhomogeneity ^[Qi] is written as a sum of three terms: first, 
Eq. (15. 6p multiplied by jJq, where we merely replace Q by Q\ and Rk by = RkTi^] 
the difference of Eqs. (15.81) and (15. 9p . both multiplied by ^, with the same substitutions; 
and Eq. (15. 10j) accompanied by |7o"rL,o- 

In order to perform the reduction of the inhomogeneity U 2 to a polynomial form, it is 
sufficient to transform the summand of Eq. (15.291) to a form involving just the polynomial 
P p (u) itself rather than its real part, i.e., half the sum of P p (u) and P p {— u). In Mellin space 
P p (u) corresponds to a polynomial in z n alone, but not (1 — z) n . As can be observed easily, 
the terms proportional to b(M) and 7o(M) are even functions of the spectral parameter 
and, as a consequence, do not change after the substitution u — > —u. On the other hand, 
the term ~ 5*2 (p) in the right-hand side of (I5.30P does not have this property. Indeed, from 
Appendix D we have 

M M 

Y,P P (-u)R P (M)S 2 (p) = J2 P p( u )Rp(M) [V 2 M - S 2 (p) - 2S_ 2 (M)] , (5.33) 

p=0 p=0 

and deduce a complimentary representation of the two-loop Baxter polynomial 

M 

Q 1 (u) = J2 p P ( u )Rp( M ) r iA M ) , (5-34) 

p=0 

with 

r ljP (M) = b{M) + 7o (M) (s x (p + M) - S 1 (M) - \ ^(p)) - \ V 2 , (p) + 5_ 2 (M) . (5.35) 

Separating the infinite-series contributions from U 2 of the form as in Eq. (15.271) . we find 
that their coefficient conspire to give the same overall coefficient but with the opposite 
sign, such that L^fQi] can be cast in the form 

Wi] 4E me[Pp - l{u)] { it Rk (M) ri ik (M) Z\ (p, k - 1) (5.36) 

p=l P ^ k=l 
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Therefore, the result for the sum of the two inhomogeneities £73 [(Jo] + ^[Qi] is free from 
non-polynomialities and reads 



M 



hJ2 ^ e[Pp : l(U)] \ [X-87o + ((7o + ) 2 + « + )+87 + r llP (M)] R p (M) (0.37) 



P =i 



V 



+ MM) \\ Z 3 {p, k - 1) - 2 (i 7o + + r llfc (M)) Z^p, k - 1) 



fc=i 



+7+ Z 2 (p, k) + 8 7o + (( 7o + ) 2 + « + ) " 16C 3 , + 2 1 +Z 2 (p) 



with 6*3,0 combining the functions 6*3,0 = 4/3 + + — ^7q a. The corresponding contri- 
bution to CJ3 is then given by 

M 

Qte*\u) = Y, 2 MPk(u)}R k (M)r^(M) , 



(5.38) 



k=0 



where we give in the same order as the corresponding inhomogeneities appear in 

Eq. (15371) 



r (j P ) (M) 7q + ( (7o+) 2 + Q+ ) ] C^) _ 1 7 + 5 3 l(A;) _ _L ^(j.) 

+^ 7 + $U*) + JE V3,x,i(k) - i 7 + V 3 , 2 (k) 
-1 7 + ((7 + ) 2 + « + ) Vs,o(*) + I C*s,o W 3 (k) - ^ 7o + W3, 2 (*) • 
It is expressed via the following nested sums 

k 



1=1 

k 



Rc,m{M) 



Z h - X (r,m - 1) 



k ^ r 

w*o=E rajRr(M) E ,„ 

r=l v ' m=l 

j r 

W*)=£ E ^(M)Z 6 (r,m), 

r=l ^ ' m=l 

k r 

V ^(k) = E ^7M E R ^{M) Z b (r } m-l). 

r=l m=l 



(5.39) 

(5.40) 
(5.41) 
(5.42) 
(5.43) 
(5.44) 
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These are related to the previously introduced sums (15.141) - (15.17P via 

W a (k) = W at0 (k), V a>b (k) = V aA0 (k), V a , b (k) = V aA0 (k), V a>b (k) = V aA0 (k). 

(5.45) 

5.2.3 Inhomogeneity U\ 

Finally, we turn to the polynomial inhomogeneity (15. 19ft of the non-polynomial part of the 
three-loop Baxter function Q\ . It can be represented as a sum 

Ui [Q { r ] ] = Uf ] [Qt P) ] + U? [Qt p) ] , (5.46) 

with 

ufW P) ] = \ {Q { 2 P \u + ^) + Q [ 2 P \u - i) - 2Q^\u)) , (5.47) 
u[\Qf p) ] = - 7o + (iu + Qt P \u + i) - iu-Qf p \u + (2M + i)Q^\u)) • (5-48) 
From this, it is immediate to find the corresponding contribution to Q% to be 

M 

Q { r m {u) = J2 *e[ft(«)] MM) r^HM), (i = 0, 1) , (5.49) 



where 



k=0 



l """ , ' m (M)=I%$(M) (b(M) + W 2 (k)) - Jl jvk , (5-50) 



r 3,k 



r 3,k 



( M )=H E {^ftW - /^(M)} , (5.51) 



and 

m 



2 \ 2' 
p=l 

(5.52) 



6 Anomalous dimensions 

Making use of the explicit solution to the Baxter equation to four loops, we can immediately 
calculate its derivatives at the argument u = ±| (see, e.g., Appendix [C]) and find the 
corresponding anomalous dimensions by means of Eq. (I2.2p expanded to the required order 



of perturbation theory. The results are 

7o = 2Si, (6.1) 



7 All results are given for even values of M. An analytical continuation to complex values can be found 
in [27]. 
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7 (asy) 



71 = -S 3 - 5_ 3 + 25_ 2 ,i - 2S 1 (5 2 + 5_ 2 ) , (6.2) 

7 2 = 5 5 + 35- 5 - 25_ 3 5 2 + 25- 2 5 3 - 245-2,1,!,! - 65_ 4 ,1 - 65_ 3 ,2 " 65_ 2 , 3 

+ 125_ 3 ,i,i + 12S_ 2 ,i,2 + 125_ 2 , 2 ,i + {S 2 + 2Sl) (35_ 3 + 5 3 - 25_ 2 ,i) 
+ (85_ 4 + 5! 2 + 45 2 5_ 2 + 25 2 2 + 35 4 - 125_ 3)1 - 105_ 2 , 2 + 16S_ 2 ,i,i) , (6.3) 

= 4 5_ 7 + 6 5 7 + 2 (5_ 3 ,i, 3 + 5_ 3 ,2,2 + £-3,3,1 + £_2,4,l) + 3 ( — 5_ 2 ,5 
+ 5-2,3,-2) + 4 (S'_2 I 1 I 4 — 5-2,-2,-2,1 — 5-2,1,2,-2 — 5-2,2,1,-2 ~~ ^l, -2,1,3 
_ 5*1,-2,2,2 — 51,-2,3, l) + 5 ( — 5_ 3 ,4 + 5-2,-2,-3) + 6 ( — 5s,_2 
+ £l,-2,4 — 5-2,-2,1,-2 — 51,-2,-2,-2) + 7 ( — 5-2,-5 + 5_ 3 ,_2,-2 
+ 5_2,-3,-2 + 5_2,-2,3) + 8 (5-4,1,2 + 5-4,2,1 — 5_5,_2 — 5_4,3 
— 5-2,1,-2,-2 + 5i, _2,1, 1,-2) + 9 53,_2,-2 — 10 5i,_2,2,-2 + H 5_3,2,-2 
+ 12 ( — 5_6,1 + 5_2,2,-3 + 5l,4,-2 +54,_2,1 +54,i,_2 — 5_3,i,i,_2 —5-2,2,-2,1 



-5i, 1,2,3 


-£1,1,3,-2 - 


£l,l,3,2 - 


-£l,2,l,3 - £l,2,2,-2 - 


£l,2,2,2 


- £l,2,3,l - £l,3,l,-2 


- 5i, 3 ,l,2 


- £1,3,2,1 - 


52,-2,1,2 


- 52,-2,2,1 - £2,1,1,3 


- 52,1,2. 


,-2 - 52,1,2,2 


~~ 52,i, 3 ,i 


— 52,2,1,-2 " 


~ 52,2,1,2 


— £2,2,2,1 — £2,3,1,1 ~ 


£3,1,1,- 


2 — 53,1,1,2 — £3,1,2,1 



,2,1,-2,1 



~ £3,2,l,l) + 13 52,-2,3 — 14 52,-2,1,-2 + 15 (52,3,-2 + 53,2,-2) 
+ 16 (5_4,i, -2 + 5_2,l,-4 — 5_2,-2,l,2 —5-2,-2,2,1 —5-2,1,-2,2 —5-2,1,1,-3 

— £l, -3,1,2 —£l,-3,2,l —£l,-2,-2,2 — £2,-2,-2,1 + £-2,1,1,-2,1 + £l,l,-2,l,-2 
+ £l,l,-2,l,2 + 5l,l,-2,2,l) — 175_5,2 + 18 ( — 54,-3 — 56,1 + £1,-3,3) 
+ 20 ( — 5i,_6 — 5i,6 — 54,3 + 5_5,i,i + 5_4,_2,1 + 5_3,_2,2 + 5-2,-4,1 
+ 5_2,-3,2 + £l,3,3 + £3,1,3 + £3,3,1 — £l,l,-2,3 — £l,2,-2,-2 — £2,1,-2,-2) 

— 21 53,4 + 22 (5i,_2,-4 + 52,2,3 + 52,3,2 + 53,_2,2 + 53,2,2) + 23 ( — 5_ 3) _ 4 

— 55,2 + 52,-2,-3) + 24 (—5_4,_3 + 5i,-4,_2 — 5i,-3,i,-2 — £1,1,1,4 — £1,1,4,1 

— £l,3,-2,l — £l,4,l,l — £3,-2,1,1 — £3,1,-2,1 — £4,1,1,1 + £-2,-2,1,1,1 + £-2,1,-2,1,1 
+ £l,-2,-2,l,l + £l,-2,l,-2,l + 5i ,1,-2,-2,1 + 5i,i,i,_2,-2 + £l, 1,2, -2,1 + £l, 
+ £ 2 ,l,l,-2,l) + 25 52,-3,-2 + 26 ( — 52,5 + 5l,4,2 + 52,4,1 + 54,1,2 + 54,2,1, 

+ 28 (£1,2,4 + £2,1,4 — £-3,1,-2,1 — £-2,1,-3,1 — £1,-2,1,-3) + 30 5_3,i,_3 
+ 32 (5i,5,i + 55,1,1 — 5_3,_2,i,i — 5_2,-3,i,i — £1,-3,-2,1 — £1,-2,-3,1 

— £2,2,-2,1 + £1,2,-2,1,1 + £2,1,-2,1,1 — £1,1,1,-2,1,1) + 36 (£1,1,5 + £1,3,- 
+ £3,1,-3 — £1,1,-3,-2 —£1,1,-2,-3 —£1,1,2,-3 — £1,2,-2,2 —£1,2,1,-3 -£2,1, 
—£2,1,1,-3) +38 5_3,_3,i + 40 (— 5i, -4,1,1 — £2,-3,1,1 + £1,1,1,-2,2) 

— 41 5 3 ,_ 4 + 42 (-5 2 ,_ 5 + 51,-4,2 + £1,-3,-3) + 44 (5i,- 5 ,i + 5 2 ,_ 3 , 2 + 5 3 ,_ 3 ,i) 
+ 46 5 2 ,2,-3 + 48 51,1,-3,1,1 + 60 (5i, 1 ,_ 5 - 5i,i,- 3 , 2 ) + 62 5 2 ,- 4 ,i + 64 5i,i,i,_ 3 ,i 
+ 68 (5i,2,-4 + 52,i,-4 — 5i,2,-3,i — £2,1,-3,1) — 72 5i,i,i,_4 — 80 5i, 1,-4,1 
-C 3 5i(53 - 5- 3 + 2 5- 2 ,i). (6.4) 



1,-3 

-2,2 
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These agree with expressions found from explicit calculations of Feynman diagrams at one 
[I] and two loops [15], and three [17] and four-loop [25] results obtained with the use of 
numerical solution of Bethe equations and the principle of maximal transcendentality. 



7 Five-loop dressing and reciprocity 

Finally, let us partially address the five loop order, namely, the one stemming from the 
dressing phase Q(u). The dressing part of the five- loop Baxter polynomial can be written 
in terms of contributions with decreasing transcendentality as 

Qf{u) = c 3 (ggV) + q[: p) ) + CsQgV) , (7.i) 

where we decomposed the term accompanying £3 according to the nomenclature of poly- 
nomial and non-polynomial inhomogeneities. Their calculation in the Wilson basis echoes 
the one performed in the previous section and yields for polynomial contributions, 

Q < fJ(u) = a^ 5 Q Q + ^S 1 T m , (7.2) 

Q%(u) = a<g 3 Q + c ( i,o) T (lj0 ) + c (2 ,o) T {2 ,o) + c (0 ,i) T m + c {l)l) T (M) + c (0)2 ) T (0 ,2) • (7.3) 

Here 

C(i,o) — — \S\ — + Si t iSi , (7.4) 
C(2,o) = |<Si — \S- 2 Si , 

c (0,l) — — f^l + — \S 2 Si — — JqSs + §Sl,-2 + gSip + 1^2,1 ) 
1 Q 2 

C(0,2) = 96 , 

and the degree-reducing constants are 

a4,c 5 = -!Si(S-2 + S 2 ) , ^ (7.5) 

ct4 ^ 3 = — Tr^i 4:SiS^ — S^S± + — -rS^Sf + S2S1 + SSi^S^ — S^^S^ + 

C G 2 o q o2 1 c C c2 o C C c2 1 G 2 c i 1 c2 c i C c2 c i 

— O3O-L — 00_ 2 OlO^ T 02010^ — ZOiOi^D^ — 2 — 2 ^ ' 2 2 1 >~> - 2 >-> 1 >~> 1 t 

— S—4S1 + -7S-2S2S1 — S4.S1 — S-2S2S1 + JS-3S-2 + 2S-3S2 + -7S—2S3 + 

ilnn lCC lOG ICC ICC ICC ICC 

+ 4>J2>J3 — _ 2^' 2 ^' 1 - 2 ~~ 2 S^ 1 ! 2 — 2^ 2 ^ 1 ' 2 ~ 2 2 2 ^ _ 2 *- 32 '-' ,2 > 1 ' 

Finally, the non-polynomial term, obeys the equation 

^B[Q { ;: l)) ] = -zS 2 1 U ( M (7-6) 

with 

tt , n 1 Q (u + i) -Q Q (u) Q (u - i) - Q (u) 

U C3 [Qo} = ~ + (7.7) 
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v 

p=l ^ k=l 

The polynomial Q** can be computed according to the method spelled out above. It 
reads 

Qt P \u) = 2Sl £) MPM\ E ^f^j E Rk{M) . (7.8) 

p=0 Z=l / fc= i 

Substituting these findings in Eq. (12.21) expanded to fifth order in the 't Hooft coupling, 
we find the dressing part of the five-loop anomalous dimensions of twist-two operators, 

7f y) = --- + C3 7 5 C3 + C 5 75 C5 , (7-9) 

with spin-dependent functions 7g 3 and ^ obeying the principle of maximal transcenden- 
tality (as well as the absence of —1 indices in nested harmonic sums) 

75 5 = §5_4 — — y5_ 3) i — 55_2,2 — §5*1,-3 + §5*1,3 + f 5 3 ,1 + 

+105'_ 2i i j i + 55i _ 2 ,i , (7.10) 

4 3 = - 5_ 6 + 5 6 + 11 5-5,1 + 4 5_ 4 ,-2 + 20 S_ 4 ,2 + 9 S_ 3 ,-3 + 14 S_3,3 + 

+4 S„ 2 ,„4 + 4 S. 2A + 3 S lt - 5 - 3 51,5 + 4 5 2 ,_ 4 - 4 5 2 , 4 - 3 5 3 , 3 - 2 5 4 ,_ 2 + 
—4 54 5 2 — 3 55,i — 32 5-4,1,1 — 8 5-3,-2,1 — 10 5_3,i _2 — 28 5_ 3 ,i,2 — 28 5_ 3 , 2 ,i + 

—6 5_2,-3,l — 2 5-2,-2,2 — 10 5-2,1,-3 — 10 5-2,1,3 — 6 5_2,2,-2 — 12 5_2,2,2 + 

-10 5_2,3,i - 20 5i,_ 4 ,i - 6 5^-3,-2 - 24 5i,_ 3 , 2 - 8 5 1 ,_ 2 ,_ 3 - 10 5i,_ 2 , 3 - 4 5 1 , 1 ,_ 4 + 
+4 5i,i,4 — 2 5i,2,- 3 + 2 5i,2, 3 + 2 5i, 3 ,_2 + 2 5i, 3 ,2 + 2 5i, 4 ,i — 14 52,- 3 ,i — 2 52,-2,-2 + 
—10 52,-2,2 — 2 52,i,- 3 + 2 52,i, 3 + 2 52, 3 ,i — 2 5 3 ,_2,i + 2 5 3 ,i,_2 + 2 5 3 ,i,2 + 2 5 3 ,2,i + 
+4 54,i,i + 36 5_ 3 ,i,i,i + 4 5_2,-2,i,i + 8 5_2, 1,-2,1 + +12 5_2,i,i,_2 + 12 5_2,i,i,2 + 
+12 5-2,i,2,i + 12 5-2,2,i,i + 36 5i,_ 3 ,i,i + 8 5i, -2,-2,1 + 8 5i,_2,i,-2 + 20 51,-2,1,2 + 
+20 5i,_ 2 ,2,i + 16 5 1 , 1 ,- 3 , 1 + 4 5 1 , 1 ,- 2 ,-2 + 12 51,1,-2,2 + 4 Si, 2 ,_2,i + 16 5 2 ,- 2 ,i,i + 
+452,1,-2,1 -2451,-2,1,1,! - 165i,i,_2,i,i. (7.11) 

7.1 Parity invariance 

As we have seen in Section 2, the conserved charge H2 acquires perturbative corrections 
which shift the bare total conformal spin of Wilson operators by their anomalous dimension 
7(g) to the renormalized one, 

j = M + l - j = M + l + lj(g). (7.12) 

This phenomenon implies that the anomalous dimensions can be defined more naturally 
as functions of the renormalized rather than bare spin, such that one can define a new 
function P of argument J, 

7 = P(M+i 7 ). (7.13) 
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The parity invariance property of anomalous dimensions is then formulated as invariance 
of P(j) under the reflection map J — > —J with J 2 = jo(jo + 1) [SHI [30]. This condition 
results in an infinite number of relations for coefficients accompanying odd powers of the 
Lorentz spin in the large-M expansion in terms of corresponding even powers. For finite 
M this property gets translated into the presence of parity-even combinations of nested 
harmonic sums flk lt k 2 ,... on \y> which can have positive even and negative odd fcj's [311 [32]. 
These functions are defined as follows. For a single-index harmonic sum, they do coincide 
with usual harmonic numbers, while for more than one index, they are defined recursively, 

— Ski > Qk lt k 2 = UJ k 1 (Qk 2 ) , ^fei,fe 2 ,fe 3 — UJ k 1 (Qk 2 , fc 3 ) ; ••• ; (7-14) 
with the involved map defined by 

UJ ki(Sk 2 M,—) = SkiMM,... ~ 2^sigii(*;i)sigii(&2)(|fci|4-|fc2]),fe3- ■ (7.15) 

Let us establish the parity invariance property for the dressing part of the five-loop 
twist-two anomalous dimensions we determined above. Expanding Eq. fl 7 . 1 3 1) in the 't 
Hooft coupling and taking into account that dressing appears firstly at four loops, we find 
the five loop dressing contribution to P(M) 

n = T5-iT4 7i-i7 4 7i, (7-16) 

where the hat on symbols denotes that we consider only their dressing parts. Explicitly, 
we use 

71 = 25-1, (7.17) 

74 = 5i(5_ 3 -5 3 -2£_ 2)1 )C3, (7.18) 

75 = C3 7 5 C3 +C5 7 5 C5 - (7-19) 
A long but straightforward calculation gives 

h ^Jt-^JV- (7-20) 

where 



P| 5 = -§ (fi_ 4 - 20i, 3 - 2fi 3 ,i - 811-2,1,1 - 4fii, 



-2,1, 



P§ 3 — — |fi_6 — 2f2_4 5 _2 — 2f2_2,-4 + 1^3,3 + 

+2 (4fi_ 4j i il + Q_ 2l l j3 + fi-2,3,1 + 4 °1 -4,1 + 2^1-2,3 + -4 + ^1,3,-2 + 

+1^3,-2,1 + 1^3,1,— 2) + 4 (0_2 ) _ 2,1,1 + 2f2_2,l ) _2,l + 3^—2,1,1,-2 + 21^1,— 2,-2,1+ 

+2Qi _2,i -2 + ^1,1,-2,-2) - 8 (3^i,-2,i,i,i + 2fii ,1,-2,1,1) + 

+2( 2 (-0-4 + (Oi )3 + fi 3il ) + 2 (317-2,1,1 + 2fi 1> _ 2>1 + ^1,1,-2)) + 

+7C4^i,i, (7.21) 

which indeed obeys the parity-invariance properties spelled out at the beginning of the 
section. 
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8 Conclusions 



In this work we have developed an improved formalism for the analytical solution of the 
multiloop Baxter equation. As a demonstration of the efficiency of the framework, we 
found the four-loop Baxter polynomial and derived in a completely analytical form the 
resulting anomalous dimensions. This was possible largely due to an improved treatment 
of superficially non-polynomial terms in the Baxter equation. While in the previous con- 
sideration, the latter yielded multiple sums with Stirling numbers of the first and second 
kind involved and, as a consequence, hampering straightforward analytical calculation of 
derivatives of the Baxter function at fixed points which enter the definition of anomalous 
dimensions. They were treated making use of the principle of maximal transcendental- 
ity [20] by writing down the most general expression in terms of nested harmonic sums 
and then fitting the multiplicative rational coefficients to numerical data. In the current 
analysis this difficulty was overcome. 

Next we found a more concise representation of the polynomial contribution to Baxter 
function by using the basis of Wilson rather than continuous Hahn polynomials. The 
complexity of these expressions was reduced roughly in half. Nevertheless it should be 
stated that our analysis of the non-polynomial terms is still in favor of a representation in 
continuous Hahn polynomials. It might still be interesting to completely restrict also the 
non-polynomial terms to a Wilson basis. As such it would be possible to compare one- 
to-one the analytic properties of the Baxter function of twist-two and -three operators, in 
order to pin down the origin of the asymptotic character of the Baxter equation as well as 
the Bethe ansatz. 

Finally, we provided further evidence towards parity invariance of twist-two multiloop 
anomalous dimensions by calculating the dressing contribution to the five-loop result and 
showing that they obey the same theorem as was established earlier. 

What we did not address are the wrapping effects in twist-two operators, which emerge 
starting from four loops. The latter are known to be described by a generalized Luscher 
formula which reads for the case of twist- L operators [2T] . 

T*~>fa> = -*ss £ a (f^) Wi + (92(1+3,) • ( "> 



It is written in terms of 



R{z, n) = Q {\z - i(n - 1)) Q (\z + |(n - 1)) Q (\z + |(n + 1)) Q {\z - |(n + 1)) , 

(8.2) 



and the function 



n-l 



l( m ~ \ n ) ~ N l( m + 1 - 5") - ¥\ 

These correspond to one-loop corrections in a sigma model. The explicit formula for twist- 
two operators was found in Ref. [21] and reads 

7 (wrap) = S 2 _ AS _ 2(3 _ ^ + ^ + 43 _ 4 «^__ 2 + AS _ 2 ^ _ 8S _^_ 2a) . (8 . 4) 
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Our method can be viewed as a first step towards the analytical computation of the 
five-loop twist-two anomalous dimensions. Wrapping contributions are not included here. 
However the knowledge of the asymptotic prediction derived with it will allow one to ana- 
lyze its analytical structure in the complex spin M plane and constrain potential wrapping 
structures. These questions will be addressed elsewhere. 
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0757394 (A.B.) and by the Russian Foundation for Basic Research through Grant No. 07- 
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A Inhomogeneities 

The inhomogeneities of the four-loop expansion of the Baxter equation appear with the 
following multiplicative functions 



Qo 



192 Qo( 



2' 



+ i Qi(f)Qo(f) i i Qo 4) (j 

X 16 r&( i \ 16 n ( i\ 192 



5 + 



e + = i 



QUI) 16 Qo(f) 192 Qo(|) 

giO) 



2> 
I I i' 



Qo(f) ' 

and analogous expressions whose argument has reverse sign being related to the above via 

or = («+)* = a+ = a, (3~ = -(p+)* = -p+ = -(3, 
X- = (X + T = X + = X, S~ = (5 + y = 5+ = 5, 

e- = -(e+)* = -e+ = -e. (A.2) 

They can be written in terms of harmonic sums making use of the explicit solution to the 
Baxter equation and read 

a = -S* + 5_ 2 , (A.3) 
(3 = —ji(Si — S3 + 35*1^2 — 35 , _2,i) , 

X — |(5-4 — S4 — 5-3,1 + 5-2,-2 — 5'_2,2 — 51,-3 + 5*1,3 — 52,-2 + 5*2,2 + ^3,1) 
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+'3-2,1,1 + "31,1,-2 — <3l,l,2 — 31,2,1 — "32,1,1 + 251,1,1,1 
+ i7o(5-3 — S3 — S-2,1 + 5'i,_2 + "Sj.,2 + <32,l) , 

8 = ai, 

e = -70(^-2 + S 2 ) + 2(5_ 3 -3 3 - 25_ 2 ,i) + 2a 1 £' 1 , 
expressed in terms of nested harmonic sums 

M 



S, 



a\,...,a k 



„, C_1 N )sign(ai) \sign(a 2 ) 

s ai ,..., ak (M) = e 1 J c e- ' 



£1=0 



£2=0 



«i 2 

t 2 



4-i 

E 

4=0 



-1) 



sign(a fe ) 



'ft 



(A.4) 



B Continuous Hahn and Wilson polynomials 

Continuous Hahn polynomials are defined as [33J 

/ , js - n (0 + c)n (a + d) n ( — n, n + a + b + c + d— 1, a + iu 
Pn (u,a,b,c,d) = z« 3 F 2 I^ a + c , a + d 

They obey the functional relation (P(u) = p n (- ■ ■ )) 

(c - iu)(d - iu) P{u + i) + (a + iu)(b + iu) P(u-i)+ (B.2) 
(2 u 2 — i (a + b — c — d) u — n 2 + (—a — b — c — d+l)n — ab — cd) P(u) = , 



1 . (B.l) 



or 



(c-iu)(d-iu) [P{u + i) - P{u)} + {a + iu){b + iu) [P{u - i) - P{u)} + (B.3) 
-(n 2 + n (a + b + c + d - 1) P(w) = . 

Wilson polynomials are defined as 

W n (u 2 , a, b, c, d) = (a + b) n (a + c) n (a + d) n 



X4P3 



-n, n+a+b+c+d— 1, a + iu, a — iu 



a + b, a + c, a + d 

They obey the functional relation {P{u) = W n (u 2 , • ■ ■ )) 

(c — iu){ia + + «)(zci + u){2u — i)P{u + z) 
-(c + m)(«-ia)(« — «6)(« — id)(2u + «)P(« — i) + 2zmAP(u) = 0. 



(B.4) 



(B.5) 



with 



A = 2 u 4 + u 2 (-An 2 -4(a + b + c + d-l)n + a + b + c + d + 
-lab - 2ac - 2bc - 2ad - 2bd - 2cd) + 
—n 2 + (—a — b — c — d + l)n — abc — abd — acd + 2abcd — bed . (B.6) 
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It can also be rewritten as 

n(n + s- l)Y{u) = B{u) [Y{u + i) - Y{u)} + B{-u) [Y{u - i) - Y(u)\ , (B.7) 

with 

n , s (u + ia)(u + ib){u + ic){u + id) , N 

B\U) = — T . (D.O) 

y ' 2u{2u + i) v 1 

The relation 

/ 1 1 1 1\ / 1 1 1 \ . . 

Pm I u, -, -, -, - I = constant x W M /2 I u, -, -, -,0 1 , (B.9) 

is immediately proved since the two recurrence relation are equal for the above choice of 
parameters. Looking at the normalization, one can check that the constant is indeed equal 
to one. 

The comparison of the two-loop Baxter equation and equations obeyed by polynomials 
yield the identifications for the continuous Hahn 

n = M , a = c=\ + ^g + \g 2 1 ^ b = d = \ - j=g + \g 2 lo , (B.10) 

and Wilson polynomials, 

n=\M , S = | + I^ 2 70) a = \, d = 0, (B.ll) 

b=1 2 + T2 9 + ^o, c=|-^ + ^ 2 7o, (B.12) 

respectively. 

C Derivatives of the Baxter functions 

The n-th derivative of Q^(u) evaluated at u = i/2 is a combination of harmonic sums 

with uniform transcendentality n and multi-indices containing only 1 and ±2. The first 
cases are 

g' (i) = -2i5 1 , (C.l) 

Q£(f) = -4 (S. 2 -S 2 + 2S ltl ), (C.2) 

gS'(i) = -24 i (S_ 2 ,i - 5 1; _ 2 + 5 li2 + S 2A - 2S 1A:1 ) , (C.3) 

Qo ^(|) = 96 (5_2,-2 — 5_2,2 — 5*2, -2 + 5*2,2 + 25-2,1,1 — 25'i i _2 i i + 25'i ) i i _2 + 

-2S lil)2 - 2^,2,! - 2^2,1,1 + AS!,!,!,!), (C.4) 

Q^i^) = — 960 i (5_2,-2,l — S-2,1,-2 + S -2,1,2 + 5-2,2,1 + 5*1,-2,-2 — $1,-2,2 + 
— $1,2,-2 + <Sl,2,2 + 5*2,-2,1 ~ 52,1,-2 + 52,1,2 + 52,2,1 _ 25_2,1,1,1 + 

25i,_ 2 ,i,i — 25'i,i,_2,i + 2jSi,i,i,_ 2 — 25*1, 1,1, 2 — 25*1, 1,2,1 — 25*1,2, 1,1 + 

-2^2,1,1,1 +45i,i,i,i,i), (C.5) 

Q*o\^) = — 5760 (5_2,-2,-2 — 5_2,-2,2 — 5_2,2,-2 + 5_2,2,2 _ 52,-2,-2 + 
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+52,-2,2 + 52,2,-2 — ^2,2,2 + 25-2,-2,1,1 — 25-2,1,-2,1 + 25-2,1,1,-2 + 
—25-2,1,1,2 — 25-2,1,2,1 — 2o_2 ) 2,l,l + 25i ; _2,-2,l — 25^-2,1,-2 + 
+201^2,1,2 + 2(31^2,2,1 + 2(31^^2,-2 ~~ 25i,i,_ 2 , 2 — 25i,i, 2 ,_2 + 25i, 1,2,2 + 
+2S'i j 2,-2,l — 2(S"i ,2,1,-2 + 25i 2,1,2 + 25^2,2,1 — 25*2^2,1,1 + 25 2 ,i,_2,l + 
— 25*2,1,1,-2 + 252,1,1,2 + +252,1,2,1 + 252,2,1,1 + 45-2,1,1,1,1 — 451,-2,1,1,1 + 

451,1,-2,1,1 — 45i,i, 1,-2,1 + 45i, 1,1, i,_2 — 45i,i,i,i,2 — 45i,i,i,2,i — 45i,i,2,i,i + 
— 45i, 2 ,i,i,i — 45 2 ,i,i,i,i + 85i, 1,1,1,1,1). (C.6) 

Let us write the two-loop Baxter polynomial as 

Qi(u) = fll (M) Q (M) + SQ 1 (u), (C.7) 

where a>i is the same as in Eq. (13.81) . The derivatives of Q^\u) are expressed in terms of 
the derivatives of Q(°\u) and those of 5Q^(u). The n-th derivative of 5Q^(u) can be 
expressed in terms of harmonic sums with uniform transcendentality n + 2. The first cases 
are 

5Qi(f) = 0, (C.8) 
«jgi(i) = -2 i (25_3 - 35_2,i - 5i,_ 2 - 5i, 2 - 5 2 ,i) , (C.9) 
8Q'{{{) = 8 (5_3,i + 5-2,2 - 5i,- 3 - 5i, 3 - 5 2 ,- 2 - 25 2 , 2 - 5 3 ,i + (CIO) 

—25_ 2 ,i,i + 25i,_ 2 ,i + 25i,i,_ 2 + 25i,i, 2 + 25i, 2 ,i + 25 2 ,i,i), 
SQ'l'(i) = -24 z (25_4,i - 5_3,- 2 + 55_ 3 , 2 + 5_ 2 ,_ 3 + 5_ 2 , 3 - 5 2 ,_ 3 + (C.ll) 
+35 2 ,3 + 53,- 2 + 353, 2 + 254,1 ~ 85-3,i,i — 45_ 2 ,_ 2 ,i + 25_ 2 ,i,- 2 + 
— 65_ 2 ,i, 2 — 65_ 2 , 2 ,i — 45i,i,3 — 25i, 2 ,_ 2 — 65i, 2 , 2 — 45i,3,i — 25 2 ,i,_ 2 + 
— 65 2 ,i, 2 — 65 2 , 2 ,i — 453,1,1 + 105_2,i,i,i — 25i,- 2 ,i,i + 25i,i,_2,i + 
+65i,i,i,_ 2 + 65i,i,i, 2 + 65i, 1,2,1 + 65i, 2 ,i,i + 65 2 ,i,i,i), 
5gS 4) (f ) = 384 (5_3,-3 - 5_3,3 - 5 3 ,_ 3 + 5 3 , 3 - 5-3,-2,1 + 25- 3) i, 2 + (C.12) 

+25_3, 2 ,i + 5_2,1,3 + 25_2,2,2 + 5-2,3,1 — 5i,_4,i — 25i,-3,2 + 

— 5i,-2,-3 + 5l,2,-3 — 25i,2,3 — 25i,3,2 — 5i,4,i — 52,-3,1 — 52,-2,2 + 
+52,1,-3 — 252,1,3 — 352,2,2 _ 252,3,1 + 53,-2,1 — 253,1,2 — 253,2,1 + 

— 35_3,i,i,i — 5_ 2 ,i,_2,l — 35-2,1,1,2 — 35-2,1,2,1 — 35-2,2,1,1 + 35i, -3,1,1 + 
+25^-2,-2,1 — 5i,_2,l,-2 + 25i,_ 2 ,i,2 + 25i, -2,2,1 + 5i, 1,-3,1 + 5i,i,- 2 ,2 + 

— 5i,i,i,_3 + 35i, 1,1,3 + 5i,i,2,-2 + 45i,i, 2 ,2 + 35i, 1,3,1 + 5l 2,1,-2 + 
+45i,2,l,2 + 45i,2,2,l + 35i,3,i,i + 52,-2,1,1 + 5 2 ,i,i,_ 2 + 452,1,1,2 + 

+452,i,2,i + 452,2,i,i + 353,i,i,i + 45_2,i,i,i,i — 45i,_ 2 ,i,i,i + 

— 45i, 1,1, i,_2 — 45i,i,i,i,2 — 45i,i,i, 2 ,i — 45i,i,2,i,i — 45i,2,i,i,i — 452,i, 1,1,1) ■ 
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D Mellin transformation 



In this appendix we devise a very efficient formalism for finding solutions to the second- 
order finite difference Baxter equation focusing on non-polynomial inhomogeneities. The 
results that we will present are very general, covering all possible powers of non-polynomial 
contributions, (u ± )~ k Qi(u), and are therefore applicable to any order of perturbation the- 
ory. The restriction and application of the machinery to three- and four-loop non- poly- 
nomial inhomogeneities, which are the main objective of the current study, are given in 
sections 15.11 and 15.21 respectively. We also employ it to the dressing-induced five-loop 
non-polynomial part in section [7J Below, we first introduce the Mellin transform for the 
Baxter function and accompanying (inverse) polynomial dressing factors in the spectral 
parameter. Then we perform a variable transformation which is particularly convenient 
for carrying out Mellin convolutions and inverse transform. 

To start with let us introduce the Mellin transform of the main ingredients. For the 
Baxter function Q(u), where the loop-order subscript is dropped being irrelevant, it takes 
the form 1281 



Q(u) =K / <koLo m -'Q{-u) 



K 



viiu r(i 



(D.l) 



While for the function with the argument shifted by i, i.e., Q(u ±i), we find that in Mellin 
space this yields a multiplication by a power of uj 



Q{u ±i) = -K du uj iu ~ l {uj t1 Q(-uj) } = K 



duuj 



iu—l 



{-uj^Q(-uj)} . (D.2) 



Next, turning to terms where Q(u) is accompanied by a positive power of a polynomial in 
the spectral parameter of the form (u + bi) L = X L , we get at first 



POO 

[u + bi) L Q(u) = X L K dujuj iX ~ 1 uj b Q(-uj) 

Jo 



(D.3) 



Then, we can re-express the product X L u lX 1 by a differential operator acting on the 
exponent, i.e., —i(d/dcouj) Ij Lj' tX ~ 1 , and subsequently integrate by parts to find 

L 



[u + bi) L Q{u) = K 



dujJ u - 1 uj- b 



IUO- 



d 

du 



{uj b Q(-u)} . 



(D.4) 



Similarly, for the Baxter function with the shifted argument Q(u±i), we get the expression 

L 



{u + U) L Q{u±i)=K 



du co iu - x u- b 







luJ 



d 
du 



(D.5) 



In the same vein we can consider terms with inverse powers of the same polynomial, 
i.e., (u + bi)~ L . First, notice that for L — 1, its Mellin integral representation 

-b 



1 



u + bi 



±i 



du uj 



iu—l—b 



0(±l=Fw) 



o 



+i j^duuj^- 1 - 
-tj^duuj^- 1 -' 



for b^O, (D.6) 
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is given in terms of the the Heaviside step function 9(x), defined conventionally as 

n/ \ f 1 , f° r X > ; / 1—. r-r\ 

y J \ , for x < . v 1 

The other observation is the well-known fact that the Mellin transform of a product of 
functions is realized as a convolution of their Mellin transforms. Namely, if a function 
Mj{u) is expressed in terms of its Mellin transform 

poo 

Mj{u) = / duuj^Fjiuj) , (D.8) 
Jo 

then the product of M 1 (u)M 2 (u) corresponds to 

M^M^u) = ^ 'dtuLU™- 1 f°° — !- Fi(cji) F 2 (—) , (D.9) 
Jo Jo ^1 

which we will write formally as 

M^M^u) M ^ r** Fl ( Ul )F 2 (^) , (D.10) 

Jo Vi \UiJ 

where the symbol — > implies the inverse Mellin transform. Putting the two results to- 
gether, leads to 

; V_ 6 , for 6^0. (D.ll) 



u + hi 



Now, upon differentiating both sides of this equation w.r.t. b, we obtain the final result 
Q(u) -> ^ ; ^ ; ' , for 6^0. 



(D.12) 

z Q(u±i) (D.i3) 



And for the shifted Baxter polynomial Q(u±i), this is replaced by 
L\ 



(u + hi) 

+ H) L+1 K j: % (ln^) L (p-y [-u^Qi-u,)] 



for 6^0. 
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D.l Change of variables 

To proceed further, we introduce the variable z as in [28J, 

"=j4l, (D.14) 
such that the Mellin transform ( ID. II) takes the form 

Q(u) = ^ dz z iu - 1/2 (1 - z)-™- 1 ^ ^) s ( D 15 ) 

where we introduced a new function (see Ref. [23J) 

Q(z) = y/z(l-z)V(z), (D.16) 
with the corresponding normalization factor 

Kl = r(i-i U )r(i + i«)- (D - 17) 

The latter leads to the unit normalization factor in Eq. (13. ip . 

Absorbing all terms of the all-order Baxter equation beyond one loop into inhomo- 
geneities Q(u) on the right-hand side, we can write it as 

{u + ) 2 Q{u + + {u~) 2 Q{u -%)- t Q {u)Q{u) = Q{u) . (D.18) 

While the corresponding equation in Mellin space reads (we omitted the overall factor 
Kiy/ z(l — z) on both left- and right-hand sides) 

z(l - z)¥'(z) + (1 - 2z)¥(z) + M(M + l)tt(z) = $(z) , (D.19) 

and where the primes on ^(z), e.g., ^f'(z), correspond to the differentiation w.r.t. the 
variable z. 

The functions Q(u) and ^(z) admit perturbative expansions 

oo oo 

Q(u) = 9* Qt(u) , *(z) = Yl ^ Mz) , (D.20) 

e=o £=o 

with their coefficients taking the form 

M M 

Q e {u) = R t,k(M) Pk(u) , * t (z) = R eA M ) z " ■ (D.21) 

k=0 k=0 

In the following, we will reserve the convention Rq^M) = Rk{M) for the leading order 
term, with the functions Rk(M) and Pk{u) introduced in Eq. (15.31) . 
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Note that the solution of (lD.19j) can be written as a sum of the z k - and (1 — z) 
expansions, i.e., 



fc=0 



i2<, fc (M) 



(D.22) 



which corresponds to the real part of its Mellin transform counterpart (1D.21j) . i.e. 



M 



Qt(u) = ^3fe[P fc (u)] R e , k (M) 



(D.23) 



fc=0 



and is the solution that is sought for. 



D.2 General properties 



Consider the polynomial solution (ID. 211) to the Baxter equation (1D.19j) . The left-hand side 
of the latter can be cast in the form 



Af-l r 

E 

k=0 



(k + l) 2 R ljk+1 (M) - (M - k)(k + 1 + M)R i)k (M) 



(D.24) 



and contains the maximal (M — l)-th power of the ^-variable compared to the solution 
(1D.2ip which is a M-th order polynomial. 

It is convenient to introduce the new coefficients r^ fc+1 (M) as 



rLk{M) 



R k (M) ' 

with Rk(M) given in Eq. ( 15. 3ft . Then, since it obeys the recurrence relation 

(k + l) 2 



(D.25) 



R k (M) = R k+1 (M) 



(M - k)(k + 1 + M) 



;D.26) 



we can replace ( ID. 241) by 



M-l 



E *"Rk{M) (k + lf 



k=0 



r ltk+1 {M) - r i>k {M) 



(D.27) 



At leading order of perturbation theory, the right-hand side ^f(z) of Eq. (ID.19I) vanishes, 
and the equation for r^ k is simply 



r! k+1 {M) - n k {M) = 0. 



'D.28) 



The solution is 



r^ k {M) = r,, (M) , R 0>k {M) = r ifi {M) ■ R k {M) 



(D.29) 
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where r^o(M) is a constant. The possible presence of nontrivial contributions from the 
(1 — z) k part of the expansion in the right-hand side of Eq. (ID. 191) does not alter any of 
the above results because its left-hand side is symmetric under the exchange z — > (1 — z). 

Note that in order to have the solution of Eq. (ID. 191) in the form of (1D.21I) . the inho- 
mogeneities Q{u) and ^(z) should be expanded in the same basis of functions 



M-l M-l 

Q(u) = Pk(u)B k+1 , = B k+iz\ (D.30) 

k=0 k=0 

and, respectively, 

M-l . M 

Q{u) = MPk{u)\ B k+1 , = - Y [ zk + (! - *)*] B k+i , (D.31) 

fe=0 fc=0 

if the (1 — z) fc -expansion constitutes a nontrivial contribution. 

Going beyond leading order of perturbation theory, the inhomogeneities in the Baxter 
equation affect the recurrence relation for the expansion coefficients ( ID.25j) as follows 



r ek+1 (M) -n k (M) = — ^ — , (D.32) 



such that the iterative solution to it reads 

k 

n, k (M) = n, (M) + ■ (D.33) 

Here, the first term in (1D.33j) corresponds to the general solution of the homogeneous 
equation and the second one is a particular solution of the inhomogeneous equation (ID. 191) . 

D.3 Polynomial inhomogeneities 



Starting from two-loop order, the right-hand side of Eq. (ID.18j) is nonzero and contains 



power-series contribution in ~ (u ± ) fc with k < 1. The terms ~ u 1 and ~ u° induce 
polynomial contributions. The bulk of them can be calculated directly in the u-space 
(see Section 4). However, it is instructive to consider their calculation also in the Mellin 
transform z-space, because, starting from 4-loop, nonpolynomial effects from lower orders 
of perturbation theory will re-emerge through these terms as well. 

For simplicity and since this is all one needs for the present analysis, we limit our 
consideration to contributions at n-th order of perturbation theory from Baxter functions 
of one loop order lower. All other cases can be treated similarly since the only difference 
will be in the coefficients accompanying u l and u°. 
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D.3.1 Contributions ~ u 1 



For the case at hand, the corresponding inhomogeneity in the right-hand side has the 
following form (see also Eq. (5.48)) 

iu + Q x (u + i)- iu~Qt (u-i) + (2M + l)Qi (u) , (D.34) 

with a function Qx(u). After the Mellin transformation to z-space (similar to Eq. (1D.15|) ) 
and cancelling the overall coefficient K\ \J z(l — z), we find 

(l-2z)& 1 (z) + 2M$ 1 (z). (D.35) 

Let now Q\{u) and ^i(z) have the following expansions (analogous to (ID. 22ft and (1D.23|) ) 



M 



M 



Qi(w) = ^3?e[P fe ( M )] Fi(k) , § 



fc=0 



fc=0 



J^L* + (!_*)* Fi(A;) . (D.36) 



Then, the corresponding particular solutions, Qi(u) and ^1(2), of the Baxter equation and 
its Mellin analogue, can be found in terms of the expansion coefficients F\(k) 



M 



M 



Q x {u) = Y,MPk(u)] F 1 (k), 



k=0 



k=0 



J2\z k + (l-z) k F l (k), (D.37) 



where the coefficients Fi(k) are determined by F±(k) by Eqs. (ID. 321) and (1D.33|) . 

It is possible to consider the expansions in z k and (1 — z) k independently. Addressing 
the z fc -expansion first, we have from (ID. 351) and (1D.36|) 



Af-1 



z k (k + l)F!(k + 1) + 2(M - fc)Fi(ife) 



fc=0 
AT-1 



z k Rk+i(M) (k + lf 



k=0 



fi(k + l) 
k+ 1 



2/i (fc) 
A; + 1 + M 



(D.38) 



where we factored out the Rk dependence and as a result introduced new coefficients 

F\{k) = R k (M)f\(k) , Fi(k) = R k {M)h{k) . (D.39) 

Next, putting Eqs. (lD\30l) . (lD~32l) and (lD~38l) together, we get for /i(A;) the following 
representation 



a(*)=/i(o)+x; 



m=l 



AH 2/ 1 (m-l) 
m m + M 



(D.40) 



Finally, when we set Qi(tt) = Qo(u), i.e. to the leading order Baxter polynomial, the 
coefficients can be found explicitly and read 



fi(k) = Sx(k) + 2[Si(k + M) — Si(M) ) . 



/i(*) = l, 

In this way we thus reproduce the 70-part of the Q\ Baxter function in (I5.30p . 



(D.41) 
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D.3.2 Contributions ~ u° 

For the u° term, the corresponding right-hand side has the following form (see also Eq. 
(5.48)) 

Q 2 {u + z) + Q 2 {u - i) - 2Q 2 {u) , (D.42) 

in terms of a function Q2(u), which we will specify later. After the Mellin transform and 
again cancelling the overall coefficient K\ z(l — z), we find for the z-space counterpart 
of the above inhomogeneity 

jjj^jS.M. (D.43) 
Here Q 2 (u) and ^> 2 {z) have the following expansions 

M M 

Q 2 (u) = J2MPk(u)\ F 2 (k) , 9 2 {z) = - [ zk + C 1 - z f} • (D.44) 

fe=i fe=i 

As it was discussed in the previous section, the expansions (ID.44j) prevent the appearance 
of negative powers of z and/or (1 — z) in the right-hand side of the Baxter equation. Then 
the corresponding particular solutions Q 2 {u) and ^ 2 (z) to the Baxter equation can be 
found in a similar fashion 

M M 

Q 2 (u) = 5>e[P fc ( M )] F 2 (k) , = - [ zk + (1 - z t] F *(k) , (D.45) 

fc=0 k=0 

where the coefficients F 2 (k) are determined by F 2 {k). 

In the following, it is convenient to separately consider the cases of symmetric and 
non-symmetric functions Q 2 (u) and ^ 2 (z) with respect to the replacements u — ► — u and 
2 — ► (1 — z), respectively. 

Symmetric case. In this case, ^ 2 (1 — z) = ^ 2 (z) and Eq. flD.43j) is equal to 

-M*) + 7T7T*2(2). (D.46) 
Z [1 z) 

Considering only the z fc -part of the expansion, we find from (1D.46|) 

Af-l 



Y,z k F 2 (k + l). (D.47) 



fc=0 

Then, introducing new coefficients via 

F 2 (k) = R k (M)f 2 (k) , F 2 (k) = R k (M)f 2 (k) , (D.48) 
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we deduce from Eqs. (1D.30I) . (1D.32|) and (ID. 470 the following representation 



/ 2 (AO = / 2 (0) + 2^ 



/2M 



m=l 



(D.49) 



If Q 2 {u) = Qo, the coefficients read 

f 2 (k) = l, f 2 (k) = 2S 2 (k), (D.50) 

and, thus, they reproduce the Qi Baxter function, when 70 = b(M) = (see (15.301) ). 

Non-symmetric case. As in the symmetric case, we can cast it in the form of the z k - 
expansion, 



M oo M M oo 

— e * k m) = E ^ E zk = H zk E zP ~ l > ( D - 

' fc=l m=0 fc=l fc=l p=fc 



51} 



where p = m + k, however, now being an infinite series rather than a finite sum. Splitting 
the series Y^Lk ^ wo P ar ts 



00 M oo 

E = E+ E 

p=fc p=fc p=M+l 



(D.52) 



we replace the right-hand side of ( ID. 511) by 



where 



E^ 1 $ 2 (p)+ E ^^(m) 

p=l p=M+l 



Hp) = T,F 2 (k). 



k=l 



Then, Eq. (ID. 531) can be represented as 



M oo M 

E ^ [Hp) - Hm)] + $ 2 (M) E = E zV ~ l [Hp) - Hm) 
p=i p=i p=i 



+ 



As a consequence, Eq. (ID. 431) reads 

, M 



1 



p=l 



= o E K _ ' + C 1 - ^ 



$ 2 (p) - $2(M) 



$ 2 (M) 
2z(l-z) 



(D.53) 



(D.54) 



$ 2 (M) 
1 - 2 
(D.55) 



(D.56) 



34 



To proceed further we have to address the nonpolynomial term in Eq. (1D.56j) . Notice 
that if we add to the left-hand side of (1D.19j) the Baxter equation obeyed by the solution 
Qi(u) with the formal condition 70 = 0, this will yield the following contribution to the 
inhomogeneity 

^)*oW, (D.57) 
where ^o(z) is the leading order solution of the Baxter equation with the property 

M 



(D.58) 



k=l 



Put together, we now have the right-hand side in the form 



M 



p=l 



J2( zP ' 1 + ( 1 - z ) P ' 1 ) Mp)-MM)+2K 2 R k {M) + 8 2 (M) + 2A 2 



1 



2z(l - z) 
(D.59) 

Choosing K 2 = — $ 2 (M)/2, we cancel the unwanted term ~ l/(z(l — z)) and find that the 
resulting Baxter equation is obeyed by the function 



Q 2 {u)-\MM)Q 1 {u)\ 10=G . 
which admits the following Mellin-space transform 



M 



p=l 



\p-i 



$ 2 (p)-$ 2 (M) (l + RpiM) 



(D.60) 



(D.6i; 



Next, using Eqs. rtHSOjl . (ITX321) and rtOTHTji we find the coefficients / 2 (Jfe) to be 

A' 



/ 2 (fc) = / 2 (0) + 



B,, 



where 



J m 2 i? m (M) 



S m = $ 2 (m) - $ 2 (M) 1 + R m (M) . 



(D.62) 



(D.63) 



Finally, taking the solution for Qi(w) from Eq. (I5.30p with 70 = 0, we find the solution for 
Q 2 (u) in the form of Eq. (ID.62I) with 



B m = $ 2 (m) - $ 2 (M) - b(M) R m (M) 



(D.64) 



D.3.3 Another form for Q\{u) 

As we explained earlier, the solution to the Baxter equation in the Mellin space is a sum 
of z k - and (1 — ,z) fc -expansions (see Eq. ( 1D.21|) ). However, the application of this form 
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as the input to inhomogeneities on the right-hand side of the Baxter equation is not very 
convenient. It is a better choice to cast it solely in terms of the z fc -expansion, 



M 



(D.65) 



fc=0 



For the four-loop calculation, considered in this study, the non-polynomial inhomo- 
geneities contain only the functions ^o(z) and ^i(z). Moreover, the function $?o(z) and 
the b(M) as well as the 70 part of ^i(z) are symmetric under the interchange z «-> (1 — z) 
and do not represent a difficulty on their own. Thus it is necessary to address the issue 
of non-symmetric contributions only for the remainder of the function ^\(z), with the 
coefficient ~ R k (M) S 2 (k). 

Our goal is to find a relation of the form 



1 M 

fc=0 



M 



z k + (l-z) k 



R k (M) S 2 (k) = J2 zk R2AM) + K 3 V ( 



(D.66) 



fc=0 



where we also added the solution of the homogeneous Baxter equation with a coefficient 
K 3 . To start with, we note that the left-hand side of Eq. (1D.66I) is a solution of the Baxter 
equation with the following right-hand side 



M 



fc=l 
M 



U* -1 + (l - «) fc_1 Rk(M)- 



(D.67) 



The second term here, i.e. ^2 k=0 (1 



\fc-l 



Rk{M), leads to the solution of the Baxter 



equation in the form ^2 k=0 (1 — z) k Rk{M) S 2 (k). Thus, we merely have to replace the 
(1 — 2;) fc -expansion by the ^-expansion. Indeed, one can show, that 



M 



At 



fc-1 



Ed 

fc=i 



— z 



k-l 



(D.68) 



fc=i 



m=l 



and Eq. (ID. 671) can be represented as 



M k 
fc=l m=l 

Thus, the general solution to Eq. (1D.67j) admits the form 

, M M 



(D.69) 



z fc + (l-z) fc R k (M)S 2 (k) = J2 zkR k(M)V 2fi (k) + K 3 ^ (z) 



(D.70) 



fc=0 



fc=0 



where V 2 fi(k) is given by Eq. f)5.40p in the main text and the coefficient K s is fixed from 
numerical agreement between the left- and right-hand sides of Eq. (1D.72j) 



-2S_ 2 (M) . 



(D.71) 
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The corresponding w-space form is then given by 

M M 

^®e[P k (u)) R k (M) S 2 (k) = ^P fc (n) R k (M) V 2)0 {k) + K 3 Q { 



u 



(D.72) 



fc=0 



k=0 



D.4 Non-polynomial inhomogeneities 



Let us now address non-polynomial inhomogeneities emerging in the Baxter equation. To 
start with, we split the Baxter function (and its Mellin transform) at £-th order as 



Qi(u + ai) = R ifi (M) + Qi(u + ai) 
with Rofl(M) = 1, where 



*t.(z) 



R e ,o{M) + ^e(z) 



M 



M 



Q(u) = Rk(M)P k (u) , = K x MM) 



z 



(D.73) 



(D.74) 



k=X k=X 

Then the inverse Mellin transform of a generic inhomogeneity is given by 



L! 



(«+) 



+ \L+X 



Qi{u + i) + (-1) 



L+X 



—l 



L+ j Qe(u-i 



(D.75) 



\ L / \ -b 
00 \ ( UJ \ 



- i L+1 K 



UJx V u i) 



-b 



-uJx Tl Q e (-uJx) 
cdx Tl Qe(-^x 



After the substitution (lD.14j) for uj —>■ z, we have for the r.h.s. of Eq. (1D.75|) 

L 



{-i) L+l K x y/z(l-z) 



— (\n Z ^ ~ Zl * 



z Jo z i 



zx(l~z) 



+ 



zi(l - z) 

(1 - z) Jz iX-zx) \"z{l-zx) 



dzx 



In 



;d.76) 



The Baxter function Q(u) and its Mellin transform *&(z) possess symmetry properties 
preserved to all order of perturbation theory: 



Q{-u) = Q{u) 



*fl -z) = V(z) . 



(D.77) 



Thus, in all calculations we will focus on just one term in Eq. ( ID. 761) . say, the first one 
on the right-hand side, and then get the complete result for the Baxter function obeying 
the required symmetry properties (1D.77j) by adding contribution with reflected argument. 
With this argument in mind, Eqs. (ID. 751) and (1D.76|) can be rewritten as 



-Q e (u + i) + (-l 



L! 



it 



-\L- 



-Q e (u-i) 



M- 



(D.78) 
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2{-i) L+1 K 1 



z(l -z)} 1/2 f z d Zl ( z(l-zi) 



. In 



To calculate the integral in the right-hand side of OP .780 . it is convenient to re-express the 
integrand involving the logarithm as a rational function of its argument via 



In 



z(l-zi) 
zi(l-z) 



d 

de 



L / z(i- Zl r 

zi(l-z) 



(D.79) 



£ = 



Now, using the expansion for ^e(z), after some algebra we obtain 



dz\ ( z(l — Z\, 



Zi \zAl-z) 



p=l L y v ^ ' k=l 



[DM) 

The first term on the right-hand side of this equation is obviously polynomial, but the 
second one is not. Note, however, that for p > M we can represent it as 



T(p + l-e) ^ * v ; k\ 



(DM) 



p=M+l yl ~ ' " ' k=l 

where the inner sum in the last term can be cast in the form 

R ^ M ) fc!r(- £ ) = Qt(u = i+ie) = Q l (u = i + is) - R , S (M) . 
Thus, we have 



(D.82) 



r( P )r(i- £ ; 



r(Jfe-e) 



fc! 



- E 



p=M+l 



T{p + l-e] 



Qi(u= § + is) 



(D.83) 



To get back to the integrand in question, we have to differentiate both sides w.r.t. e and 
set it to zero. Then one immediately finds that the last term becomes 



k=0 p=M+l F 

where we introduced notations for 



d 

de 



L-k 



r(p)r(i-e) 

r(p + 1 - e) 



;d.84) 



e=0 



L! 

fc!(L- A;)! 



Qt [u = ±f) , Qf (±|) = Q*(±|) , (D.85) 
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and for the coefficients 



de 



T(k+1 -e) 
T(p+l-e) 



(D.86) 



e=0 



which are expressed in terms of harmonic numbers as follows for a few low values of L, 

Z (p, k) = l, 

Zi(p, k) = Si(p, k) , 

Z 2 (p, k) = Slip, k) + S 2 (p, k), 

Z 3 ip, k) = Slip, k) + SSiip, k)S 2 ip, k) + 2S 3 ip, k) , 

Z 4 (p, k) = Slip, k) + 6£?(p, k)S 2 ip, k) + SS^p, k)S 3 ip, k) + 3S 2 2 (p, k) + 6S 4 (p, k) , 
Z 5 (p, k) = Slip, k) + lOSfip, k)S 2 ip, k) + 20S 2 (p, k)S 3 (p, k) + 15^(p, k)S 2 2 ip, k) 



where 



+ WSiip, k)S±(p, k) + 20S 2 (p, k)S 3 ip, k) + 2AS 5 ip, k) 
S L ip,k) = S L ip)-S L ik) 



(D.87) 
(D.88) 



and Z L (p) = Z L ip, 0). 

Assembling all results, we have 



\L+1 



' Q t (u + i) + !t^Qt(v-i) - 



( U +)L+1 



Z— ✓ pi pL 

p=l 



1 P_1 

-1V E 



(L-l)! 



E 



fc=i 



RkAM) 
k 



Z L -i(p, k — l) 



E 



fc!(L-jfe)!^ V2y — p 

fc=0 v ' p=M+l r 



Z L -kip) 



(D.89) 



And after returning to the w-space the final result, free from Stirling numbers, reads 



1 ~ (-D L+1 ~ 

Qiiu + i)+ )J +1 Qeiu-i) 



\L+1 



-2(-i) 



L+l 



M 



E 

p=l 

- E 



p 



p 



+ 



i 

- iv E 



(L-l)! 



k=l 



RkA M ) 

k 



Zl-i(p, k - 1) 



Z L _ fc (p) 



(D.90) 



^ fc!(L - fc)! - p 

where P p _i(u) is given in Eq. (15.31) . Note that for L > 0, Eq. (ID. 901) can be simplified to 
1 ~ (-1) L+1 ~ 

( M+) L+1 Q^ U + + ( M -)L+1 ^ - Z ) 
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-2(-z) 



L+l 



L 

- £ 



1 M 



p * — ' k 

p=l ^ fc=l 



1 



Z L _ k (p) 



P 



(DM) 



k=0 v ' p =M+l 

Repeating all above calculations, it is possible to obtain in an analogous fashion 



u 



\L+1 



Qe(u) 



(D.92) 
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~>2{-i) L ^K lS /z{l-z) 



1 " z^ 1 ^ fl M (Af) 
p=i ^ fc=i 



+ £ 

fc=0 



k\(L-k) 



Q 



(*), 



i) E ^^*(p) 

p=M+l ^ 



and as a consequence 
1 



2(-z) 



L+l 



M 

77 E 



"- Ll " , E /.',,( .!/) Z, (/,/,-) 



p=i ^ fe=i 



+£ 



k\(L-k)\^ e 1 ^ 2^ p 
Last but not least, as a particular case of (ID. 921) and (1D.93[) . we find 

1 ("1) L+1 M~ 



Z L ~k{p) 



p =i 



P 



1 



+ 



1 



,£+1 



(«") 



L+l 



2(-i) Z+1 E^^^(p) 



(D.93) 



2(-^ L+1 ^ v^T^) V ^(p) , (D.94) 



(D.95) 



E On the degree reducing constants 

We have seen that the polynomial-type contribution to the Baxter functions takes the form 
of a linear combination of higher derivatives T( a ,b,c,... ) plus a constant times the leading order 
Baxter polynomial Qq. This constant is fixed by the requirement that the monomial u M 
cancels in the total contribution. Thus, the following ratios are all we need to determine 
the degree reducing constant 



coefficient of u M in T< 



P(a,b,c,...) 



(a,b,c,... ) 



coefficient of u M in Qq 



(E.l) 
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Here, we list the cases which are needed for computing a 2 , and a 3 . The notation is 



S a ,K.. = S a , h ,...(M), S aA ... = 5 a , 6 ,...(2M), S aA ... = S aA ...(M/2). 

P(i,o) — 2Si - 3S±, 
P(o,i) =25*2, 

P( 2 ,o) = 95 2 - 125151 + 45 2 + 55 2 - 45 2 , 
P(i,i) = -6SxS 2 + 45 x 5 2 - 25 3 , 
P(0,2) = 125*2 + 125 4 , 

p m = -275? + 545i5 2 - 36S?Si - 45.S2.S1 + 365 2 5i + 

+85 3 - 185 3 + 305 2 5i - 245i5 2 + 165 3 , 
P( 2 ,i) = 18525? + 125 3 5i - 245 2 5x5i + 25 2 + 85 2 5 2 + 

+105 2 5 2 + 65\ - 85 3 5i - 85 2 5 2 , 
p (1>2) = -365i5 2 2 + 245i5 2 - 245 3 5 2 - 365i5 4 - 245 5 + 245 4 5i, 
p (0 ,3) = 1205 3 + 3605 4 5 2 + 2405 6 . 

In the case of 6 3 , we need the ratios 

P(0,0,2) — 252, 

P(o,o,4) = 125 2 + 125 4 , 

P(o,o,4) = 125 2 + 125 4 , 

P(o,o,4) = 125 2 + 125 4 , 

,0(0,0,6) = 1205 2 3 + 3605 4 5 2 + 2405 6 , 

P(o,2,o) =45 2 + 25 2 , 

P(o,2, 2 ) = 85 2 5 2 + 165 3 5i + 45| + 125 4 , 
P(o,2, 2 ) = 85 2 5 2 + 165 3 5i + 45 2 + 125 4 , 

p (024) = 245 2 3 + 485 2 5 2 + 1925i 5 3 5 2 + 1685 4 5 2 + 965 3 2 + 485 2 5 4 + 

+1925i5 5 + 2405 6 , 
P(o,3,o) = — 85 3 — 125 2 5i — 45 3 , 

P(o,3,2) = -165 2 5 3 - 485 3 5 2 - 245 2 2 5i - 725 4 5 x - 325 2 5 3 - 485 5 , 

P(o,4,o) = 165 4 + 485 2 5 2 + 325 3 5i + 125 2 + 125 4 , 

P(o,4,o) = 165 4 + 485 2 5 2 + 325 3 5i + 125 2 + 125 4 , 

p (0 ,4,o) = 165 4 + 485 2 5 2 + 325 3 5i + 125 2 + 125 4 , 

P(o,4, 2 ) = 325 2 5 4 + 1285 3 5 3 + 96 5 2 5 2 + 2885 4 5 2 + 2565 2 5 3 5i + 

+3845 5 5i + 245^ + 645 3 2 + 1685 2 5 4 + 2405 6 , 
p (0 5 0) = -325f - 1605 2 5 3 - 1605 3 5 2 - 1205 2 2 5i - 1205 4 5i - 805 2 5 3 - 485 5 , 
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p (060) = 645? + 4805 2 5? + 640535? + 7205 2 2 5 1 2 + 7205 4 5? + 9605 2 5 3 5i + 

+57QS 5 S l + 1205 2 3 + 1605 3 2 + 3605 2 5 4 + 2405 6 , 
P(i,o,o) — Si — 25i, 
P(i,o,2) = — 45 4 5 2 + 25 4 5 2 — 25 3 , 

P(10 4) = -245i5 2 2 + 125i5 2 2 - 245 3 5 2 - 245 4 5 4 - 245 5 + 125 4 5i, 
P(i,2,o) = —85? + 45i5 2 — 85 2 5i — 25 3 + 25 2 5i, 

P(i,2, 2 ) = -16525? - 405 3 5? + 85 2 5i5 2 - 165 2 5i - 485 4 5 4 + 165 3 5i5i + 

-165 2 5 3 - 245 5 + 45 2 5i + 12S 4 Si, 
p (1)3i0) = 165? - 8SiSf + 365 2 5? + 205 3 5 4 - 12S 2 S 1 S 1 + 65 2 2 + 65 4 - 45 3 5i, 
P(14 0) = -325? + 165! 5? - 1285 2 5? - 1125 3 5? + 485 2 5 1 5 1 2 - 725 2 5! + 

-725 4 5i + 325 3 5i5i - 405 2 5 3 - 245 5 + 125 2 5 4 + 125 4 5i, 
P( 2 ,o,o) = 45? - 45i5! + 5? + 25 2 - 5 2 , 

P(2,o,2) = 85 2 5? + 85 3 5 4 - 85 2 5i5i + 65 2 + 25 2 5? + 65 4 - 45 3 5 4 - 25 2 5 2 , 
P(2,2,o) = 165? - 165! 5? + 45?5? + 325 2 5? - 45 2 5? + 165 3 5! + 

-165 2 5 1 5 1 + 65 2 + 25 2 5? + 65 4 - 45 3 5 4 - 25 2 5 2 , 
P( 3 ,o,o) = -85? + 125 4 5? - 65?5 4 - 125 2 5 4 + 65 2 5 4 + 5? + 

-45 3 + 65 2 5 4 - 35 4 5 2 + 25 3 . 
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